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fH , Abstract 

Qh! An intriguing possibility is explored that the solar neutrino data could 

' be used as a probe to the magnetic field structure in the sun. Various cyclic 

- ^ , phenomena occurring on the surface of the sun have been accounted for by the 

XI 

( ' so-called dynamo-mechanism. According to this, a self-generating mechanism 

of solar cyclic activities gives rise to twisting toroidal magnetic fields in the 
convective zone. Although its magnitude is not known, the orientation of 
the twist is certainly opposite in northern and southern hemispheres of the 
sun. We show by numerical calculation that the solar neutrino flux, being 
sensitive to the twist, could exhibit observable seasonal oscillation, provided 
that the twist is sizable in magnitude and the neutrinos have reasonably large 
magnetic moments. This oscillation is ascribed to the fact that the earth's 
orbit is slightly inclined to the sun's equator, and that solar neutrinos pass 
through the northern (southern) hemisphere of the sun around September 
(March). We also argue that similar seasonal oscillation could be exhibited in 
the azimuthal asymmetry of recoiling electrons scattered by the solar neutrinos 
which are expected to be observed in the super-Kamiokande detector. 
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1. INTRODUCTION 



The deficit of tlie solar neutrinos observed in Homestake [1], Kamiokande II 
(KII) [2] and the gaUium experiments [3,4] has been one of the strongest stimuh 
for theoretical ideas in elementary particle physics. Among many others, the most 
interesting mechanisms which attracted a lot of attention to solve the solar neutrino 
puzzles are : (1) resonant amplification of neutrino oscillation proposed by Mikheyev 
and Smirnov and Wolfenstein (MSW) [5] , and (2) possibility of large magnetic mo- 
ments of neutrinos, which convert solar neutrinos into sterile right-handed ones [6,7]. 
These seminal ideas were followed by further elaboration thereof [8] and a few review 
articles [9] on these subjects are now available. 

These mechanisms have been expected to open new windows in neutrino physics, 
and to suggest what lies beyond the standard model. In the present paper, however, 
we would like to take a little bit different attitude of looking at the solar neutrino 
data: we will explore a possibility to make use of the solar neutrino data to investi- 
gate the magnetic structure of the sun, supposing that neutrinos are endowed with 
reasonably large magnetic moment. 

The observation of solar activities has a long history of its own. Everybody now 
agrees that phenomena occurring on the surface of the sun are outer manifestations 
of inner dynamics which controls the magnetic structure of the so-called convective 
zones of the sun. A task set upon solar physicists has been to find a coherent 
explanation of self-generating mechanisms of the cyclic phenomena of the sun. The 
cycle is such that the toroidal magnetic force lines undergo periodical change in 
orientation in the convective zones, which leads to the 11-year cycle of the sunspots 
occurrence. 

One of the important implications of the solar dynamo is that the toroidal mag- 
netic field receives twisting in the course of its time- development. Interesting enough 
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is that the orientation of the twisting is very hkely to be opposite in the northern 
and southern hemispheres of the sun. We will focus our eyes on this peculiar fact, 
and will investigate the possibility that the evidence in favor of opposite twisting in 
each hemisphere has already been encoded in the solar neutrino data. 

The study of the neutrino propagation in the presence of twisting magnetic 
fields was started by Vidal and Wudka [10] and also by Aneziris and Schechter [11]. 
Smirnov [12] later clarified the implications of the twist and pointed out a new 
mechanism of resonant neutrino oscillations. Qualitatively we can easily speculate 
the following phenomena. First of all let us recall that the earths's orbit makes 
an angle of 7.25° with the plane of the solar equator. In other words, neutrinos 
observed on the earth have passed through the northern (southern) hemisphere of 
the sun around September (March). This means that neutrinos would experience 
twisting magnetic fields in opposite directions depending on the season, i.e., effects 
of the twisting fields on the neutrino flux would exhibit semi-annual variation. 

The purpose of the present paper is to put these qualitative considerations on a 
more quantitative basis. Preliminary results have been reported in ref. [13]. We will 
show here that, with reasonably large magnetic moments, solar magnetic flelds and 
sizable twisting, the total amount of solar neutrino flux should go up and down every 
year. We will also argue that there should be seasonal change in the observation 
of the azimuthal asymmetry in the distribution of recoiling electrons scattered by 
neutrinos which are to be observed in the super-Kamiokande detectors. 

The seasonal variation to be discussed in this paper reminds us of the similar one 
discussed several years ago by Okun, Voloshin, and Vysotsky [6]. They were aware 
of the fact that there is a small equatorial "sht" in the solar magnetic held where 
the magnetic strength is very weak. They noticed that there would be a maximum 
in the left-handed neutrino flux in June and December. Perhaps it should be born 
in our mind that, while their ideas are interesting in their own right, the origin of 
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their seasonal variation differs from ours. In their case the latitude dependence of 
the magnetic field strength is the source of the variation, while in our case the twist 
is playing the crucial role. 

The present paper is organized as follows. First of all in sec. 2, we recapitulate 
the neutrino evolution equation to explain how the twisting magnetic field could 
affect the neutrino flux. In sec. 3, we describe the profile of the solar magnetic 
field, and set up a simple model for it which will be used in our numerical cal- 
culations in sec. 5. The azimuthal asymmetry of electron recoil in KII (and/or 
super-Kamiokamde ) detectors is discussed in sec. 4 on the basis of the model in 
sec. 3. Sec. 6 is devoted to discussions and summary. 



Let us begin with the evolution equation of neutrinos in the presence of twisting 
magnetic fields. The evolution equation is usually discussed in a parallel way with 
the neutral kaon system. In the presence of magnetic fields coupled to neutrino spin, 
the evolution equation is to be derived on the basis of relativistic wave equations. 
For a pedagogical derivation, readers are referred to Mannheim's paper [14], where 
both Dirac and Majorana cases are discussed. 

In our present investigation, we consider Majorana neutrinos, whose time evolu- 
tion is described by the Schrodinger equation 



where we include into our state vector ^, electron and muon neutrinos and their 
charge conjugation (denoted by C) 



2. EVOLUTION EQUATION AND TWIST 




(1) 




(2) 
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Inclusion of tau neutrinos are rather straightforward. To give a definite form of the 
Hamiltonian in (1), we take henceforth the moving direction of neutrinos the 2;-axis. 
Thereby the Hamiltonian becomes. 



H ^ 



H+ M 



\ 



(3) 



where 



H 



± 



Am^ sin^/4£; 



i^Am2sin^/4E ±V^ + Am'^ cos29/2E J 



(4) 



M 



-li{B^ + iBy) 



-KBx - iBy) 



\ 



I 



(5) 



The matter potentials perceived by electron and muon neutrinos are denoted by 
= G{2Ne — Nn)/V2 and = —GNn/V2, respectively, where Ne and A''„ are 
the number densities of electrons and neutrons, respectively. The lepton-number 
violating Majorana masses niee, irif^f^ and appear in the evolution equation only 
in the combination of Am^ = ^ (mgg — m^^)^ + ^m'^^, together with the mixing 
angle 9. 

In the course of deriving eq. (3), it has been assumed that the neutrino energy 
neutrino (transition) magnetic moment /x, magnetic field strength B, and neutrino 
masses satisfy E ^ mij{i,j = e,iJ,),E ^ iiB. The reason for the presence of the 
term 



Bx ± iBy = Bt exp{±i(l)) 



(6) 



is due to the fact that the off-diagonal elements in M come from the spin-fliping 
terms of the magnetic moment interactions. Terms that contain the 2;-component of 
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the magnetic field would appear in the off-diagonal part of H±. They are, how- 
ever, suppressed by a power iiB^/E and have been neglected in (3) in our approxima- 
tion. The phase (f) is often referred to as geometrical or topological phases, because 
it would give rise to Berry's phase [15] if it were adiabatically time-dependent. 

To have a physical insight into the phase in eq. (6), it is more sensible to 
truncate our 4x4 matrix into the following 2x2 matrix equation 



. d 

dz 



/ 










V 







V/2 —(iBt exp{i(, 

-liBTexp{-i(j)) -V/2 



(7) 



The potentials and have been replaced in the above by ±V^/2, where 

y = K + V;. = V2G{N, - N^). (8) 

This is allowed because wc can always take out common diagonal phase. Note also 
that the time derivative has been replaced by the ^-derivative in (7), since neutrinos 
travel nearly with the light velocity. The mass term Am^ has been neglected in 
(7) in order to make our argument as definite as possible. The combination of 
and {v^l)^ as a Dirac field has been considered by Zel'dovich and Konopinski and 
Mahmoud (ZKM) [16] and will be referred to as ZKM Dirac scheme. 

The effect of geometrical phase was first discussed in the hope that the existence 
of the phase might be efficient enough to solve the solar neutrino problem even 
with very small neutrino magnetic moment. Further investigation, however, showed 
that this hope was an illusion and that the adiabatic approximation used often 
in connection with Berry's phase was not applicable in this problem [11, 12, 17]. 
Although the original hope was turned down, Smirnov [12] realized that there exists 
a potentially interesting conspiracy between the geometrical phase and the resonant 
spin-fiip as well as resonant spin-fiavor transitions. 
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The best way to see the effect of the phase is, following Smirnov, to transfer to 
the new state basis 



I'eL 



exp(-i0/2) 

exp(i0/2) 



The change of the basis transforms eq. (7) into 



(\/ + 0')/2 
-hBt -{V + (t>')/2 




. d 


/ 






dz 


\ 







(9) 



(10) 



where 0' is the space derivative along the direction of the neutrino path. It is now 
obvious from eq. (10) that spatial variation of the phase amounts to either increasing 
(0' > 0) or decreasing (0' < 0) the potential. The resonant spin-flavor transition is 
modified sensitively in accordance with the sign of 0'. 



3. A SIMPLE MODEL OF TWISTING TOROIDAL 
MAGNETIC FIELDS IN THE SUN 



The solar activities and in particular the solar magnetic field profiles seem to 
be extremely complicated at first sight. Observations over many years, however, 
revealed characteristic patterns of solar activities. One of the most familiar facts 
is the cyclic occurrence of sunspots in every 11 years. Observational investigations 
show that the sunspot zones drift equatorially both in northern and southern hemi- 
spheres. The magnetic nature of sunspots has also been studied. The sunspots show 
up always in pair and the polarity of the pair is opposite. Notably, alternations in 
magnetic-field polarity of bipolar sunspot groups in every 11 years have been con- 
firmed and the oscillation period is now doubled to 22 years. The sunspots give 
us much information of the lower and middle latitude of the sun, while around the 
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poles the polarity of the magnetic fields are known to change at the maximum phase 
of the solar cycles. 

These observations were extremely strong driving forces to figure out the origin 
of the self-sustaining mechanism of the cyclic phenomena. An important ingredient 
to study the solar magnetic fields is the differential rotation of the sun. In plasma 
fluid with high electric conductivities, magnetic force lines tend to go along the 
motion of the fluid. The differential rotation elucidates therefore how to generate 
toroidal magnetic flelds from poloidal ones. Although the differential rotation is 
crucial to explain some of the solar phenomena, another mechanism converting the 
orientation of the toroidal magnetic flelds had to be looked for elsewhere to explain 
the cyclic phenomena. 

The global convection current has been called for to flU in the missing link of 
the cycle. The global convection is a global-scale fluid motion in the convective 
zone, extending over the sphere. Qualitatively, the Coriolis force is responsible for 
twisting the magnetic flelds and the time-evolution processes of magnetic flelds are 
illustrated in Fig. 1 schematically following the work of Yoshimura [18]. One can see 
in Fig. 1 that the toroidal magnetic flelds are twisted in the course of development. 
The reversed toroidal magnetic flelds locate deep inside the convective zone. The 
toroidal flelds in the outer area arc pushed toward the surface of the sun, spread and 
create sunspots, and then dissipate. Eventually they are replaced by the reversed 
toroidal flelds, that is the polarity reversal at the maximum phase of the solar cycle. 

A quantitative analysis of the solar cycle requires detailed study of the magneto- 
hydrodynamial (MHD) equations, which describe the time evolution of magnetic 
flelds under given velocity flelds of plasma media. To handle the MHD equations 
analytically, however, is nearly impossible and numerical integrations are the best 
means to solve the MHD equations. 

Yoshimura [18] has been engaged in elaborating numerical solutions to the MHD 
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equations. He prescribed the velocity fields of differential rotation and global con- 
vections by introducing several adjustable parameters. His numerical computations 
reproduce the basic characteristics of the solar activities. It has been observed that 
magnetic field tubes (torus) encircling the rotational axis of the sun appear natu- 
rally. The reversal of orientation of toroidal magnetic fields takes place inside the 
convection zone. The reversed magnetic fields propagate towards the surface of the 
sun, which turns out to be solar activities of the sunspots. The parameter depen- 
dence of his simulation was also examined. Although the magnetic fields around the 
circling tube is twisting, the twisting parameter is not uniquely determined. 

In the present paper, being motivated by Yoshimura's computer simulation, we 
will use a simplified model of the solar magnetic fields. We consider toroidal mag- 
netic fields winding along the tori which sit in the northern and southern convective 
areas (Fig. 2). Each torus is parametrized by its radius a, the latitude of its cen- 
ter A, and the distance b measured from the center of the sun. The latitude A in 
Fig. 2 shows the direction of the path of neutrinos which are assumed to have been 
produced at the center of the sun. 

We take the magnetic fields to be distributed symmetrically around the rotating 
axis of the sun and the twisting proceeds at the constant rate along the torus. We 
introduce a parameter X which parametrizes the degree of the twist and is defined 
as the distance along the torus to wind it once. (See Fig. 3.) The orientation of 
the twist varies in the first and second 11 years and is opposite in the northern and 
southern tori. The patterns are schematically shown in Fig. 4. 

The direction of neutrino propagation is taken to be 2;-axis and the longitudinal 
(latitudinal) direction along the x (y)-axis. In this coordinate system the twisting 
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magnetic fields are parametrized by 



By 



^ rix cos 9{z) ^ 



= B 



Uy sin^(2;) cosQ;(;i;) 
riz sin ^(2;) sinQ;(2;) 



Here angles 6{z) and a{z) shown in Fig. 5 are defined by 



.(.) = tan-M^; 



(11) 



(12) 



cosaU 



6cos(A — X) — z 



(13) 



with 



R{z) = A/[6sin(A - A)]2 + [6cos(A - X) - z] 



(14) 



The parameters Ux, riy and Uz in eq. (11) are sign factors specifying the twist which 
depends on seasons and either the first or the second periods of 11 years. They are 
given from Fig. 4 by 



(+, +, — ) for Nl: northern hemisphere, first 11 yr, 

(— , +, +) for SI: southern hemisphere, first 11 yr, 

(— , — , +) for Nil: northern hemisphere, second 11 yr, 

(+, — , — ) for SII: southern hemisphere, second 11 yr. 



(15) 



So much for the parametrization of the magnetic field. We are now in a position 
to evaluate Bt and the phase from eq. (11), and we find 



Bt = B^J cos^ 9{z) + sin^ e{z) cos^ a{z), 



(16) 
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(I){z) = ntan~^[— (6cos(A - A) - z)] + S. 



(17) 



Here a sign factor n and S are given by {n, S) — (+1, 0), (—1, tt), (+1, tt), (—1, 0) for 
NI, SI, Nil, and SII, respectively. We are interested in the z-derivative of in the 
evolution equation 



-1 
1 



2-/-'^ foJ (18) 



l + (27r/X)'l(,cos(A-A)-2]' SI, SII. 



This shows clearly that the effect of the phase depends crucially on whether neutrinos 
pass through the northern or southern hemispheres of the sun. 

Qualitatively, if neutrinos go through the northern hemisphere (around Septem- 
ber), then the negative sign on the right hand side of eq. (18) indicates suppression 
of the the diagonal element of the Hamiltonian in (3), resulting in an relative en- 
hancement of the off-diagonal effect. We can expect large spin-flip conversion rate 
around September and considerable neutrino flux suppression. On the contrary if 
neutrinos are passing through southern hemisphere (around March), then the neu- 
trino flux will be less suppressed. The seasonal contrast depends on the magnitude 
of the magnetic field. If the magnetic field strength is larger (e.g., in the maximum 
phase of the solar cycle), then the seasonal variation is rather veiled. This is qual- 
itatively different from what Okun, Voloshin and Vysotsky [6] have considered. In 
their case, larger the magnetic field, more conspicuous the seasonal variation. 



4. AZIMUTHAL ASYMMETRY OF THE ELECTRON RECOIL 



Let us turn to another possibility to detect the twist of the solar magnetic field. 
A conceivable direct evidence of the twist could be found in the azimuthal asym- 
metry in the distribution of the electrons scattered by solar neutrinos. The origin 
of the asymmetry is the interference between electromagnetic scatterings due to 
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the neutrino magnetic moment and weak processes. Such an asymmetry has been 
discussed by Barbieri and Fiorentini [19] in their search for an evidence of hehcity 
conversion. They considered Dirac neutrinos with magnetic moment. 



Neutrinos produced at the center of the sun are purely in the left-handed state, 
whose probability amplitude is given by (z/gL, (^'/^L)'^) = (1)0). In the course of 
traversing the solar convective zones, the amplitudes acquire complex phases. Sup- 
pose that the neutrinos just entering the super-Kamiokande detector are described 
by the amplitudes 

(z>eL,(i>ML)'') = (l«L|e^^Mafl|e^^«). (19) 
The spin of the neutrinos then has an expectation value 



^ I aLan \ cos{(pL - ^Pr) ^ 



< S >= 



V 



I ttLttR I sin((/7L - ipr) 
(I a« P - I |2)/2. ] 



(20) 



The phases </7l, and of the left and right components are calculable numerically 
by use of (10). We immediately notice that the particular direction of the neutrino 
spin (20) may produce an asymmetry with reference to this direction. In fact, an 
asymmetry can be found if we look at the azimuthal distribution of the recoiling 
electrons. The phases </7l and ipR are sensitive to the twisting magnetic field in the 
sun, we might be able to expect a semi-annual variation in the azimuthal asymmetry 
data. 

The differential cross sections relevant to our case becomes as follows: 

(7rr r/rr*^^ drr dn"^ drr^^ drr'^'^ 



dTdif dTdip dTdip dTdip dTdip dTdip 

Here, T is the electron kinetic energy and ip is the azimuthal angle of the recoiling 
electron. The weak gauge-boson exchange contributions for eUeL ^-UeL (e(t'^L)'^ 
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e(t'/iL)'") and electromagnetic one are denoted by (ia^/ {da'^y ) and rfo-gm, respectively. 
These cross sections together with their interference ( d(T^^^ and d(T^^^) are explicitly 
given by 



d^w_ 
dTdifi 



QlOr- 



rrieT 



El 



(22) 



dap. 



int 



dTdip 



u ,0 , 1 1 , 



dTdip hb 2mg T Ei, 



(23) 



(24) 



d(7w 

dTdif 



2 C^nie 



(25) 



1^ = (€t • k)(^)-^^b. + .i(l - 1^)]. 
dl dip hb 2v27rme-/ -C/;. 



(26) 



The couplings are listed below. 



1 1 

gL = - + sin^ 6li^, gr^ = -- + sin^ 6'^', g^ij = sin^ Ow 



(27) 



The momentum of the recoiling electron is denoted by k. Since we have neglected 

neutrino masses, these cross sections are common for Majorana and ZKM neutrinos. 
Note that the cross sections (22), (23) and (24) were considered by Barbieri and 
Fiorentini [19] and eqs. (25) and (26) are our new addition. 

We have introduced a notation the transverse component of the neutrino 
spin < S > (times 2), namely, 



^ I aiaR I cos((/?L - (pn) ^ 
I aLOR I sin((/7L - (/?r) 




(28) 
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The existence of in eq. (24) and (26) shows clearly that the azimuthal distri- 
bution of kicked electrons would exhibit asymmetry with respect to the direction 
(28). Barbieri and Fiorentini [19] have calculated (/j^ — </?i as a function of /iB for 
untwisting magnetic field case. They pointed out that the asymmetry measurement 
appears to be within the reach of observation in future experiments. The point we 
would like to stress here is that the phases Lpi and (fR are potentially sensitive to 
(f)' as we see from the evolution equation (10) and we anticipate that the axis of the 
asymmetry varies from season to season. 

Before concluding this section, we present a formula for the suppression rate of 
neutrino flux observed in KIT experiment 



/ dE^^E,; SSM) /^-- dTdif 



da 



RiKll) = — ' ^ ' , ' -. (29) 

!dE,^{E,-SSM)!^--dTd^^{\aL\^^l) ^ ' 

Here ^{Ei,; SSM) denotes the neutrino energy spectrum calculated in the standard 
solar model (SSM). The integration over T is taken from the threshold energy W of 
Kll detector up to maximally possible value T^ax = T^^[Ey). 

The suppression rate differs from | p, since {v^jll)^ also contributes to the 
detection [20]. If the magnetic moment is on the order of (10~^° ~ 10~^^)/xs, the 
electromagnetic contribution is negligibly small and this becomes 

R{Kll) =1 aL 1^ +0.144 \aR\'^ . (30) 

In passing note that the suppression rate in the Homestake experiment is simply 
given by R{C\) —\ ql since their detector is blind to (i^/iL)*". 

5. NUMERICAL ANALYSIS OF THE SOLAR NEUTRINO FLUX 
AND THE AZIMUTHAL ASYMMETRY 
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We have undertaken numerical analysis of eq. (10). The integration is performed 
from z — Zq to z — Zi, where Zq and Zi are the entrance and exit coordinates of the 
magnetic zones 

zo,i = b cos(A - A) ± ^a2-62 sm\A-X). (31) 

The initial condition is set as {fygL, (i^^uL)*") = (1, 0) at 2; = zo and various parameters 
are taken as follows: 

A = ±7°, A = 15°, a = 0.1694i?Q, 6 = 0.88137?©. (32) 

The values A, a, and b were read off from the Yoshimura's computer simulation. 
The choice of the latitude of the neutrino's path A = 7° (—7°) corresponds to the 
period from August to October (from February to April). In these latitude, the 
magnetic field strength varies only a little, and we have assumed in our calculations 
that B in eq. (11) is spatially constant. We take the matter potential [21] 

V(z) = O.178i?0^ exp[10.84^1 - z/R^jcm-^ (33) 

in the convective zone. 

One might perhaps wonder how it could be that a + 6 is greater than the solar 
radius Rq as in (32). The reason for our choice of the values a and b in (32) is as 
follows. Although we said that magnetic field tubes were formed as in Fig. 2, the 
shape of the tubes is not an exact torus according to the Yoshimura's simulation. 
Since we are thinking of neutrinos passing through in the latitude A = ±7°, we 
have to take a model reproducing the magnetic fields at the edge of the tubes in 
lower latitudes. The choice of the values a and b are taken to realize these realistic 
magnetic fields along the neutrino path. 
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We are basically interested in the fiB- dependence of the flux of the left-handed 
electron neutrinos P{i'l) —\ cll P= 1— | clr where and are the probability 
amplitudes in (19). For a better understanding of numerical computation, we note 
that at 2; = O.85ii!0, for example, 



V{z = Q.85Rq) = 1.70 X 10-^"cm-\ (34) 

= 0.85«,) = ±(1.44 X 10-")^^^^J|f||L_cm-. (35) 

This is to be compared with 

Mi^-2.9xl0-'(^)(^)cm- (36) 

In other words, if the twist parameter X is a fraction of Rq and fiB is of the order 
of (1 ~ 10) X 10~^°/iBkG, then the seasonal variation is likely to be seen in P{i^l)- 

In Fig. 6, P{i^l) is plotted as a function of fiB, for various twist, X = (tt/S)/?©, 
X — 'kRqi X = 27: Rq, and X = SttRq. The solid (broken) curves correspond to the 
cases where neutrinos go through the 7° latitude in the northern (southern) hemi- 
sphere, that is , around September (March). The solid curve denoted by "toroidal" 
corresponds to the X — oo case. 

Characteristic features discussed quahtatively can be seen in Fig. 6. When (l)'{z) 
is negative as in September, then the off-diagonal element of the Hamiltonian (3) 
is more cfi^cctive and the spin-flip process occurs more frequently than the opposite 
case. In fact the neutrino flux in September is smaller than that without twist {X — 
oo). When (f)'{z) is positive as in March, things become opposite, i.e., the suppression 
of the flux is weakened. The seasonal differences are distinct for relatively weak 
magnetic fleld = (2 ~ 6) x lO^^^/iekG. For stronger magnetic flelds, say, 
liB = (8 ~ 10) X 10~^°yLiskG, the seasonal difference is obscured by the frequent 
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spin-flip processes and even opposite tendency arises in comparison with the weak 
magnetic field. 

To confront these numerical computations with long-term observational data at 
Homestake and KII, we have to vary the magnetic field strength B in accordance 
with the phase of the solar cycle. In the light of our poor knowledge on such 
long-term variation of B, we take here a phenomenological approach. Namely, we 
evaluate B by assuming that it is proportional to i/sunspots number. In Fig. 7, 
the sunspots number in the last two decades are displayed together with our fit. At 
the maximum phase of the solar cycle, we input a value (/xSjmax , and the 11-year 
variation of B is determined by making use of Fig. 7. 

In Fig. 8, we compare our estimation of the seasonal variation with the ^^Cl data 
in the last twenty or so years. Our numerical calculation was done for two cases 

(a) X = (7r/3)i?0, ii^B)^^ = 10 x lO-^VskG, 

(b) X^ttRq, ifiB)^^ = 9 X 10-1 VskG. 

Broken lines connecting predicted values are drawn only for a guide. Our calculations 
show the seasonal variations more clearly in the quiet period of the sun than in the 
active period. Although the observational data are afflicted with large error, it is 
not unfair to say that our predicted fluxes reproduce the data rather well. One can 
see an up-down structure of the flux in the CI data. Our claim is that this structure 
could be interpreted as a manifestation of the twisting magnetic fleld. 

The KII data are also compared with our numerical calculations of (29) in Fig. 9. 
Our numerical estimate shows again an up-down structure in particular in the quiet 
period of the sun (i.e., 1986-1988). The available KII data are those in the maximum 
phase of the sun and the expected structure is not so conspicuous. It is also to be 
mentioned that the data are averaged over 7 to 9 months and the up-down structure. 
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if it were present, may have been washed away by the averaging procedure. It is 
more desirable to take data in a larger detector (like super Kamiokande) without 
taking average over long term. 

The difference of the azimuthal angle (fiR — (pL is evaluated in Fig. 10 for the 
above two cases, (a) and (b). Black (white) triangles are those around September 
(March) for the case (a) and black (white) circles around September (March) for 
(b). It is obviously seen that the asymmetry axis swings between two completely 
different directions. As far as the asymmetry measurement is possible in a large 
volume detector as was remarked by Barbieri and Fiorentini, the measurement of 
the seasonal variation in so large angles should also be feasible. It is therefore a good 
proposal in future experiments to measure the seasonal variation of the asymmetry 
axis. 

Finally we add a little comment on the adiabaticity. In our present calcula- 
tions wc have solved (10) numerically and there is no reason to worry about the 
adiabaticity. It is, however, more instructive to check the validity of the adiabatic 
approximation by using our explicit parametrization (18) and (33). The adiabaticity 
means that the mixing of states proceeds very slowly during the passage of neutrinos 
in a typical precession length. More explicitly, the adiabaticity is ensured (near the 
resonance point) provided that 

2(//S)2»l/'(^) + (/)"(z). (37) 

Numerically the parametrization (33) tells us V'{z) = (2 ~ 3) x 10~^°cm~^ for 
0.85 < z/Rq < 0.90. On the other hand eq. (18) gives (f)"{z) = (-4 ~ 4)xl0-2°cm-2 
for 0.85 < z/Rq < 0.90 and X — {n/3)RQ. These values should be compared with 
(liBf = 8.41 X 10-22(^/io-iVB)^(S/lkG)2cm-2. One can see therefore that for 
the parameter range of /iB adopted in our numerical calculations, the adiabaticity 
condition is not well satisfied. This is of course consistent with the observation made 
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previously in literatures [10, 12, 17]. 



6. SUMMARY AND DISCUSSION 

In the present paper we have been seeking for a possibility of making use of 
the solar neutrino data to probe the solar magnetic structure. More specifically 
we have examined the effect of twisting magnetic field to the solar neutrino flux 
and the azimuthal asymmetry which could be measured in the super-Kamiokande 
experiment. We have pointed out that within a reach of observation there should be 
semi-annual variation in the solar neutrino flux and the azimuthal asymmetry axis, 
provided that the twist parameter is sizable and /iB = (1 ~ 10) x 10~^°//BkG. 

There have been a lot of attempts of putting upper bounds on neutrino mag- 
netic moments from both terrestrial and celestial observations. The Dirac-type as 
well as Majorana-type magnetic moments are constrained by antineutrino-electron 
scattering data [22] (| |< 4 x 10^^°/xb), and stellar energy loss through neutrino 
pair emission [23] (| /i |< 1 x 10"^^//^). There is also a possibihty of deriving a 
bound from a reactor experiment [24] . The synthesis of ^He in the big bang gives 
us a bound for a Dirac- type magnetic moment [25]. There are several studies of 
restricting the magnetic moment more severely from supernova 1987A [26]. This 
is, however, still under debate [27]. In any way, the Majorana-type moment is not 
constrained by supernova 1987A. Considering the ambiguity of the solar magnetic 
fleld strength, we believe that our choice of the parameter /iB is within a reasonable 
range. 

ACKNOWLEDGEMENT 

The authors would hke to thank M. Ogura for collaboration in the early stage of 
the present work. 

19 



REFERENCES 



[I] R. Davis, Jr., in Neutrino 88, Proceedings of the XIII International Confer- 

ence on Neutrino and Astrophysics, Boston, Massachusetts, 1988, edited by J. 
Schneps et al. (World Scientific Pub., Singapore, 1989) p. 518. 

[2] K. Hirata et al, Phys. Rev. Lett. 63, 16 (1989); ibidem 65, 1297; 1301 (1990). 

[3] P. Anselmann et al, Phys. Lett. B285, 390 (1992). 

[4] A.L Abazov et al, Phys. Rev. Lett. 67, 3332 (1991). 

[5] L. Wolfenstein, Phys. Rev. D17, 2369 (1978); D20, 2634, (1979); S. P. 
Mikheyev and A.Yu. Smirnov, Nuovo Cimento C9, 17 (1986); Yad. Fiz. 
42, (1985) 1441 [Sov. J. Nucl. Phys. 42, 913 (1985)]. 

[6] L.B. Okun, Yad. Fiz. 44, 847 (1986) [Sov. J. Nucl. Phys. 44, 546 (1986)] ; 
L.B. Okun, M.B. Voloshin and M.l. Vysotsky, ibidem 91, 754 (1986) [44, 440 
(1986)]; Zh. Eksp. Teor. Fiz. 91, 754 (1986) [Sov. Phys. JETP 64, 446 
(1986)]. 

[7] A. Cisneros, Astrophys. Space Sci. 10, 87 (1981); K. Fujikawa and R. Shrock, 
Phys. Rev. Lett. 45, 963 (1980). 

[8] C.S. Lim and W. Marciano, Phys. Rev. D37, 1368 (1988); H. Minakata and H. 
Nunokawa, Phys. Rev. Lett. 63, 121 (1989); E.Kh. Akhmedov, Phys. Lett. 
B213, 64 (1988). 

[9] S.M. Bilenky and S.T. Petcov, Rev. Mod. Phys. 59, 671 (1987); T.K. Kuo, 
Rev. Mod. Phys. 61, 937 (1989); J. Pluido, Phys. Reports 211, 167 (1992). 

[10] J. Vidal and J. Wudka, Phys. Lett. B249, 473 (1990). 

[II] C. Aneziris and J. Schechter, Int. Journ. Mod. Phys. A6, 2375 (1991); Phys. 
Rev. D45 (1992) 1053. 

20 



[12] A. Yu. Smirnov, Phys. Lett. B260, 161 (1991); JETP Lett. 53, 291 (1991); 
E.Kh. Akhmedov, P.L Krastev and A.Yu. Smirnov, Z. Phys. C 52, 701 

(1991) ; A.Yu. Smirnov, in Proceedings of Joint International Lepton-Photon 
Symposium and Europhysics Conference on High Energy Physics (Geneva, 
Switzerland, July 1991) p. 648. 

[13] T. Kubota, T. Kurimoto, M. Ogura and E. Takasugi, Phys. Lett. B292, 195 

(1992) . 

[14] P.D. Mannheim, Phys. Rev. D37, 1935 (1988). 

[15] M.V. Berry, Proc. R. Soc. London A392, 45 (1984). 

[16] Ya. B. Zel'dovich, Dok. Akad. Nauk. SSSR 86, 505 (1952); E.J. Konopinski 
and H. Mahmound, Phy. Rev. 92, 1045 (1953). 

[17] S. Toshev, Phys. Lett. B271, 179 (1991). 

[18] H. Yoshimura, Astrophys. Journ. 178, 863 (1972); Astrophys. Journ. Suppl. 
467, 29: iMdem 52, 363 (1983). 

[19] R. Barbieri and G. Fiorentini, Nucl. Phys. B304, 909 (1988). 

[20] C.S. Lim, M. Mori, Y. Oyama and A. Suzuki, Phys. Lett. B243, 389 (1990); 
L Rothstein, Phys. Rev. D45, 2583 (1992); H.Minakata and H. Nunokawa, 
Phys. Rev. D45, 3316 (1992). 

[21] H. Spruit, Slolar Phys. 34, 277 (1974). 

[22] W.J. Marciano and Z. Parsa, Annu. Rev. Nucl. Part. Sci. 36, 171 (1986). 
[23] J. Bernstein, M. Ruderman and G. Feinberg, Phys. Rev. 132, 1227 (1963). 
[24] P. Vogel and J. Engel, Phys. Rev. D39, 3378 (1989). 



21 



[25] J. A. Morgan, Phys. Lett. 102B, 247 (1982); M. Fukugita and S. Yazaki, Phys. 
Rev. D36, 3817 (1987). 

[26] J.M. Lattimer and J. Coopcrstcin, Phys. Rev. Lett. 61, 23 (1988); R. Barbieri 
and R.N. Mohapatra, Phys. Rev. Lett. 61, 27 (1988); S. Nussinov and Y. 
Rephaeh, Phys. Rev. D36, 2278 (1988). 

[27] M.B. Voloshin, Phys. Lett. B209, 360 (1988). 



22 



Figure Captions 

Fig. 1: Schematic pictures of the twist of toroidal magnetic fields in the convective 
zone. The twist is generated by the differential rotations of the sun and the 
global convection of the plasma fluid. These pictures are given from the work 
of Yoshimura [18]. 

Fig. 2: A model of twisting toroidal magnetic fields. For simplicity, we assume that 
magnetic fields are winding along the tori which are located in the northern 
and southern hemispheres. The parameter A represents the latitude of the 
path of neutrino which is assumed to be created at the center of the sun. The 
latitudes of neutrino paths are restricted by | A |< 7.25°. 

Fig. 3: Parametrization of the twist by X which is the distance along the torus to 
wind it once. 

Fig. 4: Twists and the direction of toroidal magnetic fields. They change symmet- 
rically, depending on the hemisphere and the period. 

Fig. 5: Description of toroidal magnetic fields. 

Fig. 6: The fiB dependence of the left handed neutrino flux for various magnitude 
of the twist parameter X — (7r/3)i?o, X — ttRq, X — 27: Rq, X — SttRq. 
Solid lines represent cases where neutrinos pass through the northern hemi- 
sphere, while broken lines correspond to cases where neutrinos pass through 
the southern hemisphere. The curve denoted by "toroidal" corresponds to the 
X — oo case. 

Fig. 7: Fit of sunspot numbers which we used to obtain the time variation of the 
magnitude of magnetic fields. 

Fig. 8: Comparison between our predictions and the Cl-Ar data. Broken lines are 
drawn as a guide. The up-down structure of our predictions is the inevitable 

23 



consequences of the twist of toroidal magnetic fields. We have analyzed two 
cases: (a) X = {'k/^)Rq, {fiB)^^ = 10 x IQ-^VfikG and (b) X = ttRq, 
il^B)^^ = 9 X lO-^VskG. 

Fig. 9: Comparison between our predictions and the KII data. Our predictions 
show a structure, but the KII data do not. This is because their data are 
those averaged for 7 to 9 months which are too long to see the structure. Also 
their observation are almost in the maximum phase of the solar magnetic fields 
where seasonal difference is minimum. 

Fig. 10: Theoretical prediction of the difference of the azimuthal angle (pR — (fL 
in the last two decades. Black (white) triangles are those around September 
(March) for (a) X = {n/3)RQ, (iiB)^^ = 10 x lO-^VfikG. and black (white) 
circles those around September (March) for (b) X — ttRq, {/jiB)^^ — 9 x 
lO-^VskG. 
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